Effect of the exchange hole on the Gutzwiller approximation in one dimension 
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The Gutzwiller approximate solution to the Gutzwiller wavefunction yields exact results for the 
Gutzwiller wavefunction in the infinite dimensional limit. Implicit in the Gutzwiller approximation 
is an approximate local form of the fermion exchange hole. This approximate form is the same for 
all dimensions but is incorrect except in infinite dimensions. We implement the correct form for the 
exchange hole into the Gutzwiller approximation. We perform calculations on the one-dimensional 
Hubbard model at half-filling. They indicate that the implementation of the exchange hole already 
brings the Gutzwiller approximation into very close quantitative agreement with the results of the 
full Gutzwiller wavefunction. Metallicity as well as anti-ferromagnetism are recovered. 

PACS numbers: 71.10.Fd,71.30+h,71.10.Ca 
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I. INTRODUCTION 

The Hubbard model [H, 0, 0, 0] and its descendants 
have contributed greatly to our understanding of strongly 
correlated systems [1 Jo, 0| and in particular of the metal- 
insulator transition ^ (MIT) exhibited by them. Early 
attempts 0, H to explain the MIT were based on the use 
of a projected wavefunction due to Gutzwiller (GWF). 
The GWF is a variational method whose starting point 
is a non-interacting wavefunction, in which double oc- 
cupations and consequently charge fluctuations are sup- 
pressed. 

Exact solutions to the GWF are known in one [1, 
and infinite dimensions [l^, [ll|, • In one dimension the 
exact solution of GWF is metallic, a conclusion which 
was shown [isj to be general for finite dimensions. The 
exact solution to the Hubbard model at half-filling in 
one-dimension is insulating for all finite values of 
the interaction strength. Extended versions of the GWF 
with charge fluctuations can only account for insulat- 
ing behavior when correlations between doubly occu- 
pied sites and empty sites are incorporated (bound exci- 
tons) [15, l6, Vt\. 

The GWF is often treated via an approximation also 
due to Gutzwiller 0, 0, i, [11 (GA). The GA pre- 
dicts a MIT [1, [l^ between a paramagnetic metal and 
a paramagnetic insulator (Brinkman-Rice transition) at 
half-filling, in contradiction with the exact Gutzwiller 
solution. The GA also does not account for anti- 
ferromagnetic correlations properly whereas exact diag- 
onalizations have shown that the GWF, inspite of be- 
ing based only on projecting out double occupations, 
reproduces anti-ferromagnetic correlations remarkably 
well [11 . On the other hand the GA corresponds to 
the exact solution of the GWF when the number of di- 
mensions is infinite [13, [HI, [13, IMI. An improved GA 
has previously been constructed by Metzner [22j where 
it is shown that self-energy corrections can restore metal- 
licity. Another important study relevant here is that of 
van Dongen et al. 23] in which it is shown that metallic- 
ity can be recovered based on dimensional scaling argu- 
ments, however finite orders of perturbation theory are 



not sufficient to remove the Brinkman-Rice MIT. In two 
dimensions metallicity can be recovered 0, HI via a di- 
agrammatic summation method in which the error terms 
are estimated with high accuracy [11 . This has also been 
demonstrated numerically (27j . 

Interestingly, the exchange hole (pair-correlation func- 
tion of particles with parallel spins, defined as g{rij) = 
{nifij) where the average is over the Fermi sea) for non- 
interacting electrons at half filling, applied in GA is inde- 
pendent of the physical dimension. It corresponds to the 
exact GWF result only in infinite dimensions (See Fig. 
[IJ . The crucial point here is that in finite dimensions the 
exchange hole extends over several lattice sites, while in 
GA it is restricted to the on-site term. The central mo- 
tivation of the present paper is to show that the failure 
of GA is due to the over-simplified approximation of the 
exchange hole. 

The GA consists of taking exchange into account in 
a combinatorial fashion. The configurations considered 
obey the Pauli principle in the sense that no two particles 
of the same spin can be found on the same site. On the 
other hand no other correlation effect exists between like 
spins in the GA. Hence the exchange hole is only local 
corresponding to the infinite dimensional case (Fig. [T|). 
It is known, however (see for example Refs. [J and (28| l 
that the exchange hole in a finite number of dimensions 
has a nontrivial functional form (Fig. [T]). 

The GA as well as extensions of it enjoy widespread 
use in a variety of strongly correlated problems. The 
Brinkman-Rice transition has been used by VoUhardt 
to describe the solid-liquid phase transition in He"^ ^M- 
More recent applications of the GA include extension 
to the time-dependent case [11, implementation for the 
multi-band case (sol . [sij , ensembles with varying particle 
number (Bardeen-Cooper-Schrieffer wavefunction) [sl . 
and the calculation of matrix elements between ground 
and excited states 'ss'] . Variants of the approximate solu- 
tion have also been applied B 4 . 35, 36] in the resonating 
valence bond method [13, [H, [39|, which is based on a 
completely projected Gutzwiller wavefunction. 

In this paper we augment the GA with an improved 
treatment of the exchange hole. Our calculation sam- 
ples the exact distribution in the occupation number rep- 



2 



II. METHOD 
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FIG. 1: Pair distribution functions of systems of non- 
interacting fermions on a lattice for different dimensions at 
lialf-filling. Tlie Gutzwiller approximation uses the pair dis- 
tribution function of the inifinite limit for all dimensions. 



resentation. The hopping term, being the expectation 
value of an operator which is not diagonal in occupa- 
tion number representation is approximated in a similar 
fashion as in the original GA. Implementation consists of 
Monte Carlo sampling first used to calculate the exact 
GWF by Yokoyama and Shiba. f40| It may also be pos- 
sible to approximate the distribution with the exchange 
hole as was done for continuous systems of interacting 
fermions. [4ll, For lattices a two-site approxima- 

tion which takes into account the exchange hole has been 
proposed by Razafimandimby [4^ . 



To motivate the central idea of our approach the di- 
mensional dependence of the exchange hole at half-filling 
is presented in Fig. [TJ Due to exchange there is an effec- 
tive repulsion between nearest-neighboring (r = 1) parti- 
cles which is strongest in the case of one dimension, and 
it decreases with the number of dimensions (see also Ref. 
[4^). Further away from the origin (r > 1) the functions 
exhibit gradually decreasing oscillations. 



Via comparison with the full GWF we show that such 
a procedure brings the GA into excellent agreement with 
the full GWF. The interaction energy is exact as a func- 
tion of the variational parameter, and the hopping energy 
also shows close correspondence where full GWF results 
are available. The anti-ferromagnetic order parameter is 
also in excellent agreement between our approach and the 
full GWF calculation, while the original GA significantly 
underestimates antiferromagnetism. 



A. The Hubbard Hamiltonian and the Gutzwiller 
wavefunction 

The Hubbard Hamiltonian [l], 0, H, 01 is given by 

N N 

H = -tY^ cLcjT + UJ2^^1^il- (1) 

{ij)<7 i=l 

We will assume a system with L lattice sites and with 
Ni and particles with spins up and down respectively. 
For future reference we also define the anti-ferromagnetic 
order parameter as 




(2) 



where Sz{i) denotes the z-component of the spin at site 
i, Ci is either 1 or —1 depending on which sublattice site 
i belongs to. 

The variational wavefunction to the Hubbard model 
with which we are concerned is the Gutzwiller wavefunc- 
tion (GWF) 



l^") = exp 



-7 



(3) 



where \FS) indicates a Fermi sea of non- interacting 
fermions, the sum in the exponential counts the number 
of doubly occupied sites, and 7 is a variational parame- 
ter. For a homogeneous system the Fermi sea is formed 
by filling in the plane wave states with the lowest hopping 
energies. 



\FS) 



-'lli-'-'-lN^i 



|0). 



(4) 



Eq. ([¥]) can be rewritten in terms of sums over configu- 
rations in real-space as Q 



\FS) = L 



{g} {h} 



o«kigi 



^ikis,^ -gi _ 'SNt 



(5) 



oill-hi 



^iIn^ hi 



^glT---^gNrT'^hii---ChN^i|0)' 



where g and h denote the configurations of particles with 
spin up and down respectively. These configurations are 
such that at most only one particle of each spin can oc- 
cupy a particular site. To save space we introduce the 
notation D [k; g] for the determinants. The normaliza- 
tion of the GWF can be written 

(*|*) ^ L-(^T+^.)^^|J)[k;g]|2 (6) 

{g} {h} 

x|D[l;h]|2exp[-27i?(g,h)], 
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where D{g, h) denotes the number of double occupations 
for the particular configuration of up-spin and down-spin 
particles g and h respectively. One can define the prob- 
ability distribution 

Pgwf{sM = |D[k;g]nS)[l;h]|2exp[-27i?(g,h)]. (7) 

The exchange hole is obtained via tracing out all but 
two variables corresponding to particles with parallel spin 
in the distribution Pgwf with 7 = 0. With the help 
of Pgwf one can write expectation values diagonal in 
occupation representation. For example, the expectation 
value of the number of double occupations can be written 
as 



nit nil 



E{g,h}-PG14/F(g,h)D(g,h) 

E{g,h}^GWJ^(g,h) 



(8) 



The hopping term is not diagonal in the coordinate 
representation. For example, if a hopping of an up-spin 
particle between particular sites i and j is considered we 
have 

(*|cI^c,t|*) =i-(^T+^^'^' E®*[k;g'..T] X (9) 

{g} {h} 

I)[k; g]|S)[l; h]pexp[-7{i?(g',^.^, h) + D{g, h)}]. 

The hopping changes the configuration from g to g'ijp 
hence the determinant as well as the projection term is in 
general altered. In Eq. © g'^^ denotes configurations 
with site i occupied and site j unoccupied and g config- 
urations with site i unoccupied and site j occupied. The 
prime on the sum over the up-spin configurations indi- 
cates that only these types of configurations enter the 
summation, that allow for the hopping. 

The hopping can leave the number of double occupa- 
tions unchanged, increase or decrease it by one. Hence 
we can rewrite Eq. ^ as 

(vl/|ct^c,T|*) =L-(^T+^.)^' E®*[k;g',,.,] x(10) 

{g} {h} 

D[k; g]|S[l; h]|2exph7{2i^(g, h) + Ai^(g',^.^, g; h)}], 

where AZ?(g'jj|, g; h) denotes the change in double occu- 
pation when an up-spin particle hops from site j to i (in 
other words the configuration changes from g to g'y|)- 
We now define the estimators for hopping from site j to 
i as 

x'SwFisM = -<S)*[k;g',^.^]/D*[k;g] x (11) 

exp[--7Ai:)(g'ypg;h)], 

xg{^^(g,h) = -tS)*[l;h',^.J/I)*[l;h] X 

exp[-7Ai^(g;h'„t,h)]. 

and use these definitions to write 

Xgwf(s, h) if site j has a particle 
with spin a and 
site i does not. 







otherwise. 



Eqs. ([TT|) and ([12]) allow us to write the concatenated 
form for the estimator for the kinetic energy as 



(13) 



and write the expectation value of the total energy as 



Egwf — 



E{g,h} PawFis, h){TGWFis, h) + UD{g, h)} 
E{g,h} ^GlVi^ (g,h) 

(14) 

where the summations are now unrestricted. A Monte 
Carlo procedure can be constructed [13] to sample the 
distributions Pgwf and evaluate the expectation values 
defined in (Eq. (P])). 

The expectation value of the energy for the GWF was 
solved exactly in one dimension by Metzner and VoU- 
hardt @ . The Fermi step in this case is finite for all finite 
values of the interaction strength, hence the Gutzwiller 
wavefunction was shown to be metallic in one dimension. 
Millis and Coppersmith have later generalized this con- 
clusion to any system of finite dimensions [13]. 



B. The Gutzwiller approximation 

In the Gutzwiller approximation the determinant fac- 
tors of the probability distribution Pgwf and of the es- 
timator Tgwf (Eq- (Uni)) is replaced by configurational 
averages obtained from the non-interacting system. Our 
description of how this is done is based on Reference 0] . 

Considering only the up-spin channel one can write the 
normalization of the Fermi sea as 



r(F^|F^)T=L-^TE|D[k;g]p = l, 



(15) 



since the wavefunctions that enter are normalized 
planewaves themselves. As the sum in Eq. (|15|) is over 
all configurations of up-spin particles on the lattice, such 
that at most one particle occupies a particular site we 
can approximate each term by its average as 



|2)[k;g]p«(|D[k;g]p) 



c 



(16) 



where C^^ denotes the number of ways N-^ particles can 
be placed on L lattice sites. The down-spin particles can 
be handled similarly. This approximation results in a 
simplified probability distribution as compared to Pgwf 



PGA(g,h) =exp[-27i?(g,h)]. 



(17) 



allowing the rewriting of averages for quantities diagonal 
in the occupation number representation. For example 
the average number of double occupations in the GA can 
be written as 



(12) 



E 



ni^Uii ) = 



E{gM-PGA(g,h)i^(g,h) 
E{g,h}^GA(g,h) 



(18) 
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Approximating the kinetic energy is complicated by 
the fact that the hopping is not diagonal in the occu- 
pation number representation. Here a configurational 
average is needed for products of two determinants 
S)*[k; g']D[k; g] for particular hopping terms. Consid- 
ering only the up-spin channel one can write 

X = ^{FS\cl^c,^\FSh (19) 

= £-^T^'S)*[k;g',^.^]S)[k;g]. 
s 

One can also evaluate the average hopping over the Fermi 
sea explicitly as, 



1 



k 



V exp[ik- (Ri-Rj))], 



(20) 



where Ri and Rj denote the pair of lattice sites involved 
in the hopping, and the asterisk indicates that the sum 
be performed over occupied states only. The configura- 
tions entering the sum in Eq. (|19p are the ones with one 
up-spin electron on site j and site i unoccupied. Of such 
configurations there are C^~^^. Hence the approxima- 
tion 



S)*[k;g' ]D[k;g] « (D*[k;g',,.]I)[k;g]) =X 



L-2 
Nt-l 



c 



(21) 

can be introduced. Using this approximation the average 
hopping of an up-spin particle from site j to site i can be 
written 



c 



L-2 
N-,-1 



(22) 



Elg} E{h} Pga{^. h)exp[-7Ai^(g^^■p g; h)] 

E{g}E{h}^GA(g,h) 

Eq. (|22p is an approximate expression for the hop- 
ping energy. It allows to formulate an estimator for the 
hopping which is diagonal in the occupation number rep- 
resentation. We can then write the estimator for the 
hopping from site j to site i (the effective operator which 
is averaged in order to calculate the hopping energy) as 



Xg^(g;h) = XTexp[-7AD(g%.pg;h)], 



(23) 



X^i(g;h) = Txexp[-7Ai?(g;h%i,h)], 



where 



C 



L-2 ■ 



(24) 



We can use these definitions to write 

Xga(Sj h) if site j has a particle 



'GA 



(g,h) 



with spin tr and 
site i does not. (25) 

otherwise. 



Eqs. ([23)) and ((25)) allow us to write the concatenated 
form for estimator of the kinetic energy again as an un- 
restricted sum 



rGA(g,h)^^r^(g,h), 



(26) 



resulting in the expectation value of the total energy as 



Ega 



S{g.h} Pga{^, h){TG^(g, h) + UD{^, h)} 

E{g,h}-PGA(g,h) 



(27) 

Note that for the non-interacting system {U = 0) the en- 
ergy is exact by construction (one could also define the 
constants To- from this condition). In summary the esti- 
mator for the hopping is the product of a scaling factor 
and a factor which accounts for the change in the number 
of double occupations caused by the hopping itself. 

In summary the GA can be considered a two-step ap- 
proximation: the exact distribution Pgwf, which is com- 
plicated by the determinant factors, is replaced by the 
simpler Pga, and the estimator for the hopping Tqwf is 
replaced by Tqa- In the Tgwf the quotient of determi- 
nants, which originates from the fact that the hopping is 
not a diagonal operator in the occupation number repre- 
sentation is replaced by an average value to arrive at the 
approximation Tqa- The approximate estimator Tqa is 
diagonal in the occupation number representation. 

In the case of the GA the average energy (Eq. (|?7)) ) 
and other relevant observables can be evaluated analyt- 
ically [1, i, 0, i, [il- The resulting MIT, known as 
the Brinkman-Rice transition Q , is characterized by the 
vanishing of the expectation value of the double occupa- 
tions as well as that of the hopping energy. The latter 
can be shown to be a result of the closure of the Fermi 
step 



C. Implementing the exchange hole 

In the following we investigate the effect of the ex- 
change hole on the Gutzwiller approximation. To this 
end we substitute Pgwf for Pqa in Eq. (|27)) . We also 
scale the constants To- so that the hopping energy re- 
mains exact in the non- interacting limit. We refer to this 
approximation scheme as the GA-X. For the expression 
of the energy we can write 



EnA- 



GA-X 



S{s,h} Pgwf{^, \^){Tga-x{s, h) + UD{q, h)} 
E{g,h}^GWJ=^(g,h) 

(28) 

The estimator for the hopping Tqa-x has the same form 
as Tqa, the only difference is the scaling factor to satisfy 
the condition in the non- interacting limit. 

The expression for the average number of double oc- 
cupations as a function of the variational parameter 7 in 
our scheme is the same as in the exact GWF case (Eq. 
([5])). This does not necessarily mean that we obtain the 
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exact expectation value of this operator, however. Since 
the hopping is approximated as a function of the param- 
eters t and U, the minimization in the parameter 7 does 
not guarantee that the exact 7 is obtained. 
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FIG. 2; Absolute value of the total energy per particle calcu- 
lated in the Gutzwiller approximation: Monte Carlo results 
for various system sizes and analytical results for the ther- 
modynamic limit. The inset shows the absolute value of the 
hopping energy. 



^2- 




FIG. 3: Variational parameter 7 as a function of the interac- 
tion parameter U. 



D. Monte Carlo sampling 

To sample the distributions described in the previous 
subsection we apply the Monte Carlo method on a lattice 
due to Yokoyama and Shiba [lO] . Our system consists of 
a fixed number of sites L and up-spin and down-spin par- 
ticles iV| and respectively. Our Monte Carlo method 
consists of generating configurations with at most one 
particle of each spin on each lattice site. We generate 




FIG. 4: Hopping energy calculated via the Monte Carlo 
sampling 40] of the Gutzwiller wavefunction, the Gutzwiller 
approximation with the exchange hole, and the standard 
Gutzwiller approximation for a system of 60 sites as a function 
of the variational parameter 7. 



the configurations using two types of moves. We attempt 
moves of particles of a particular spin to sites without 
particles of that spin. We also attempt exchange moves 
between sites occupied by particles of opposite spin, and 
between sites which are doubly occupied and empty. We 
checked our Monte Carlo code against a full GWF cal- 
culation based on exact diagonalization for 12 sites and 
found excellent agreement. 



III. RESULTS 

In the following we apply the method described above 
to the one-dimensional Hubbard model at half-filling. In 
order to test our Monte Carlo program we perform cal- 
culations for GA. We have performed calculations with 
sizes up to 180. The result in the thermodynamic limit 
is also known Q. AH of these calculations are based on 
MC runs of on the order of 10^ steps. The hopping and 
interaction energies were calculated for an J7 = 1 sys- 
tem. In order to stabilize the search for the minimum, 
for all of the following calculations, we fitted the calcu- 
lated data pointsji^. The results for the absolute value 
of the total energy are shown in Fig. [2] as well as the ab- 
solute value of the kinetic energy shown in the inset. The 
Brinkman-Rice transition is clearly visible for the curve 
corresponding to the thermodynamic limit. The calcula- 
tions of systems with different sizes converge to the curve 
in the thermodynamic limit, in particular they converge 
to zero in the region of U where the Brinkman-Rice tran- 
sition predicts insulating behavior (where the energy is 
zero). 

We note that often the Fermi step is taken as the indi- 
cator of the Brinkman-Rice MIT 0, S^] . To calculate the 
Fermi step the density of states in momentum space is 
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FIG. 5: Absolute value of the hopping energy calculated from 
the GA-X approximation for three different system sizes. 



calculated. One can argue that as the Fermi step closes 
the hopping energy becomes zero, and for an open Fermi 
step the hopping energy has to be finite Q- Hence we 
can take the hopping energy as an indicator of the MIT. 
Moreover, the interaction energy, or rather the number 
of double occupations, is zero as well. 
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FIG. 6: Comparison of the antiferromagnetic order parameter 
(defined in Eq. ^) for the exact Gutzwiller wavefunction, 
GA-X, and GA. 

In Fig. [3] the variational parameter 7 is shown as a 
function of the interaction parameter U for the full GWF 
and the GA-X calculations. Quantitative agreement is 
found between the two curves. In Fig. |4] the hopping 
energies as a function of the variational parameter 7 are 
compared from an exact calculation for the full GWF, the 
GA-X, and the GA schemes for a system at half-filling 
with 60 lattice sites. In the case of 60 sites Monte Carlo 
sampling [40| was used in all cases. The agreement be- 
tween the exact GWF and GA-X is excellent, the exact 
curve essentially coincides with the GA-X results. GWF 



and GA-X differ only in the definition of the estimator 
for the kinetic energy. The fact that the kinetic energies 
of the two approaches essentially coincide is indirect ev- 
idence that the GA-X approximation is metallic. More 
evidence for this conclusion is provided by comparing the 
hopping energy for different system sizes, shown in Fig. 
[S] for sizes 36, 48, and 60. The hopping energy, whose 
becoming zero indicates the Brinkman-Rice MIT, shows 
negligible size-independence and does not become zero in 
the range of U considered (whereas the size-dependence 
is strong for the GA (Fig. 

Since the interaction energy as a function of the vari- 
ational parameter is identical for the GWF and GA-X, 
our numerical evidence suggests that the GWF and the 
GA-X coincide. It is of interest to note that the GA 
and the GWF coincide in the case of infinite dimensions, 
where the form of the exchange hole assumed in the GA 
coincide with the exact exchange hole. Our results indi- 
cate that implementing the exact exchange hole in one 
dimension brings the GA into agreement with the GWF 
in one dimension. 




FIG. 7; Absolute value of the total energy per particle calcu- 
lated via the GA-X: Monte Carlo results for various system 
sizes and the results for the thermodynamic limit of the full 
GWF (Ref. d, The inset (logarithmic in energy) shows 

that the finite size dependence of the energy increases with 
the interaction strength U. 



In Fig. [5] we present the anti-ferromagnetic order pa- 
rameter (defined in Eq. ([2])) for a system of 60 lattice 
sites. The three different methods GWF, GA-X, and GA 
are compared. Implementation of the exchange hole re- 
covers anti-ferromagnetism entirely, a results that can be 
anticipated from the results on the hopping energy (Fig. 
S]). We stress that metallicity and anti-ferromagnetism 
are recovered together as the exchange hole is imple- 
mented. In Refs. [4^ \^ expressions are given for 
the spin-correlation functions. Using these expressions 
we have calculated the antiferromagnetic order parame- 
ter defined in Eq. ^ and we have obtained a value of 
(M^) = 0.0201 for a system of 60 sites with 7 ooin 
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IV. CONCLUSIONS 
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FIG. 8: Double occupation (shown on a logarithmic scale) 
calculated via GA-X and the thermodynamic limit (Ref. 

M)- 

good agreement with our result of (Af^) = 0.022 ±0.005. 
Our recovery of metallic behaviour is related to the re- 
covery of anti-ferromagnetic behavior. 

In Figs. [7] and [H] we present a comparison of the to- 
tal energies and the double occupations between various 
system sizes. In the inset in Fig. [71 which shows the 
absolute value of the energy on a logarithmic scale, the 
energies appear to converge to the exact solution of the 
GWF in the thermodynamic limit. The logarithmic de- 
viations appear to be larger at large U . The average 
number of double occupations as a function of the vari- 
ational parameter 7 are shown in Fig. [51 The exact 
expression for the number of double occupations in the 
thermodynamic limit as a function of 7 for the GWF 
wavefunction is given in Refs. din. The agreement 
between the exact result and the GA-X results is goo for 
lower values of 7, discrepancies appear only at large val- 
ues of 7, noticible mainly at values larger than the value 
that minimizes the energy aX U = 20 (see Fig. ([3|)). 
The discrepancy can be partly attributed to finite size 
effects, which tend to oversetimate ordering, and thereby 
suppress double occupations, and also the difficulty in 
sampling with the standard MC method in systems 
which are approaching a critical point (7 00). 



We have investigated the effect of implementing the 
exchange hole in the Gutzwiller approximation in one 
dimension for the Hubbard model at half-filling. The es- 
timator used for the hopping energy was taken from the 
Gutzwiller approximate solution, but the distribution of 
configurations was the one corresponding to the exact so- 
lution of the Gutzwiller wavefunction, with the exchange 
hole implemented. Comparison with exact calculations 
based on the Gutzwiller wavefunction were presented. 

The resulting approximation is in excellent quanti- 
tative agreement with the exact result. The approxi- 
mate hopping and the hopping of the exact solution of 
the GWF are in excellent agreement. When compared 
to the Gutzwiller approximation we find that through 
implementing the exchange hole metallicity and anti- 
ferromagnetism of the full GWF are both recovered. Our 
essential conclusion here is that by implementing the ex- 
change hole one can account for the anti-ferromagnetic 
correlations present in the exact GWF, and this proce- 
dure results in the recovery of metallicity. 

As is well known pH [l3 . [2l| the Gutzwiller approx- 
imation is exact in the limit of infinite dimensionality. 
What our results suggest is that implementation of the 
exchange hole brings the Gutzwiller approximation into 
agreement with the full GWF results. In the future we 
plan to study this question in more than one finite di- 
mensions. 

From a methodological point of view our result may 
lead to improvements in the future. It is possible 
to construct approximate potentials for the exchange 
hole [U [43, [U . Such an approximate potential could 
then be sampled which would result in an approximate 
GWF, or the exact distribution could be sampled facil- 
itated by the approximated potential via umbrella sam- 
pling [ill, stochastic potential switching [ioj or acceler- 
ated Monte Carlo methodsfs^,!!!!!!,!!!. 
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